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Abstract 

This paper deals with correction of shape distortion due to zero-mean normally dis- 
tributed errors in structural sizes which are random variables. A bound on the maximum 
improvement in the expected value of the root-mean-square shape error is obtained. The 
shape correction associated with the optimal actuators is also characterized. An actuator 
effectiveness index is developed and shown to be helpful in screening actuator locations in 
the structure. The results are specialized to a simple form for truss structures composed 
of nominally identical members. The bound and effectiveness index are tested on a 55-m 
radiometer antenna truss structure. It is found that previously obtained results for optimum 
actuators had a performance close to the bound obtained here. Furthermore, the actuators 
associated with the optimum design arc shown to have high effectiveness indices. Since only 
a small fraction of truss elements tend to have high effectiveness indices, the proposed screen- 
ing procedure can greatly reduce the number of truss members that need to be considered 
as actuator sites. 


Introduction. 


Manufacturing errors are recognized as a major cause of shape distortion in large space 
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structures (e.g. Ref. 1). These errors, as well as many other sources of distortion (such as 
variations in the coefficient of thermal expansion), are random in nature, and a statistical 
analysis needs to be carried out to estimate the resulting shape distortion (e.g. Refs. 1-3). 
When actuators are used to reduce shape distortion, the corrected shape is also random in 
nature arid requires a more complex and costly statistical analysis (Ref. 3). 

Trusses are one candidate for large space structures. When these trusses are used as 
backup structures for electromagnetic reflectors, maintaining accurate surface shape is im- 
portant. Optimizing locations of actuators used to reduce shape distortion on large trusses is 
a formidable problem because of the discrete nature of the problem. Actuators are typically 
located in truss elements, and given in truss elements and n actuators we have (”) possi- 
bilities of locating them. Several ad hoc methods for tackling this combinatorial problem 
were proposed (e.g. Refs 4-6). However, these methods are not guaranteed to converge even 
to a local optimum, and therefore the quality of the actuator configuration they produce is 
difficult to evaluate. 

The present paper has two objectives. The first objective is to show that under some 
reasonable conditions it is possible to obtain a lower bound to the expected value of the 
root-mean-square (rms) of the corrected distortion. The second objective is to develop an 
effectiveness index for truss members that can be used to screen them as potential sites 
for actuators. This is accomplished by characterizing the correction shape associated with 
optimal actuators, and comparing the correction shape associated with an actuator located in 
a given truss member with the optimum shape. Results obtained by a heuristic optimization 
algorithm in Ref. 3 are used to validate the bound and effectiveness index. 

Linear Distortion Correction. 

The structure is assumed to have a set of n a sensors that measure n a components of the 
distortion field, and a set of n actuators used to correct the distortion. Denoting the vector 
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of sensor-measured components of the distortion as 0, we assume that the actuators seek to 
minimize a weighted root- mean-square (rms) measure of the distortion 0 mis defined as 

V’rms = (1) 

where B is a positive semideiinite weighting matrix. 

It is convenient to work in a transformed (modal) coordinate system where the matrix 
B is the unit matrix. That is we write 0 as a linear combination of modes 0, i = l,...,n s 

(A = $0 (2) 

where the matrix $ is composed of columns 0,- which are orthonormal with respect to B. 
That is 

$ r Z?3> = I (3) 

where I denotes the unit matrix. Then from Eqs. (1) - (3) 

0L» = 0 T 0 (4) 

We assume that the behavior of the structure and actuators is linear, so that if a unit action 
of the zth actuator produces a displacement vector tq, then the corrected shape vector 8 is 
given as 

n 

8 = 4' + ^ tq Qi = 0 + UO (5) 

i=i 

where 0 is a vector of actuator outputs (with components 0 t ) and U is a matrix with columns 
U{. Note that the components of all of these vectors arc modal (amplitude) coordinates with 
the physical vector (denoted with a hat) obtained by the modal transformation, c.g., 

u, = <K (6) 

The rms of the corrected distortion is 

SLs = PBS = 8 t 8 ( 7 ) 
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The optimum vector 0 winch minimizes 6 rilis is easily shown (e.g., Ref. 4) to be the solution 
of the system 

AO = r (8) 

where 

A = U T U (9 a) 

r = -U T x/> (9b) 

and then the corrected shape 6 is 

6 = 4>-UA~'U t 4> (10) 

and 

■ 5 .™ = ^,,. -VUA-'V* ( 11 ) 


Statistical Analysis. 

Let the distortion field be due to a set of errors or disturbances in the structure char- 
acterized by their amplitude parameters e, r = Often the statistical properties of 

the c,’s are known and we want to obtain the statistics of <p and 6. Since the behavior of the 
structure is linear the total distortion due to the N errors, V’> is given as 

A 

0 = ^ f.0. = Ve (12) 

i=i 

where ti\ is the shape distortion due to a unit c,-, 'I' is a matrix with i/>, as its ?th column, 
and c is the vector of e,’s. 

We assume that the e^’s arc zero-mean, normally distributed random variables with a 
covariance matrix E. From Ecp (12) the covariance matrix of 0 is 

Cm = FW T ) = yE(ee T )<H = (13) 
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where E denotes the expected value. Using Eq. (4) we get that the expected value of ip 2 m3 


is 




WL.) = 


(14) 


t=l 


that is, the trace of Cw The expected value of 8 2 m3 is obtained by rewriting Eq. (11) as 


8 2 = ib 2 - W T W 

u rmj Vrms rr 


(15) 


where 


W = L T U T xl> 


(16) 


and L is a square factor of A -1 (i.e., A 1 = LL T , e.g., the Cholesky factor). The covariance 
matrix of W is then given as 

(17) 


Cmv = L t U t C^UL 


and 




- E^nv),, 


(18) 


1 = 1 


The effectiveness of the actuators is measured by the distortion ratio, g, defined as 


# 2 = ^(C.) 


(19) 


&«,.) 

Optimum Actuators 

The problem of optimizing the properties of the actuators (such as location) requires the 
minimization of E(6 2 ms ) which is the same as the maximization of the trace of Cww- We 
rewrite Eq. (17) as 

Cww — U’ C^U" ( 20 ) 


where 


U m = UL 


( 21 ) 


First we note that U' consists of a set of orthonormal vectors. Indeed 


{U’ T U')~ l = ( I/U t UL ) _1 = ( L t AL )~ 1 = L~ l A~'L~ T = L~ l LL T L~ T = I (22) 
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It is shown in Ref. 7 (theorem 1) that where Xj are the eigenvalues of C ^ arranged as 

n, 

max ^(CW),-,- = ^ Xj (23) 

t j=n s -n+l 

where Xj are the eigenvalues of C^, arranged as 

0 < A, < A 2 . . . < A na (24) 

We are able to realize this maximum by taking the columns of U" to be the eigenvectors 
Vj, j = n s — n + 1, ...,n corresponding to these largest eigenvalues. + From Eq. (21) we see 
that this is equivalent to taking U as a set of linearly independent combinations of these 
eigenvectors. Altogether we have that 

E(6L.U = - E A > < 25 > 

1=1 j = n 5 ~n + 1 

For most problems it will be impossible to select actuators that produce displacement fields 
which are exact linear combinations of the eigenvectors Vj corresponding to the largest eigen- 
values of C^. However, it may be possible to come close. In particular, we can judge the 
suitability of an actuator by checking how close is its displacement field to a linear combi- 
nation. If the eigenvectors v 3 of C^, are normalized to unit length 

vjv 3 = l (26) 

then an effectiveness index of the ith actuator (which produces the displacement field u,) is 
*7 i 

<7.=( ( 27 ) 

j=n$ — w + 1 

which is the cosine of the angle between u, and the hyperplane spanned by the last n 
eigenvectors. Only actuators with an effectiveness index close to one are likely to be good 
actuators. 

+ From the derivation in [7] it is also clear that any linear combination of these eigenvectors 
will do. 
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Applications to Truss Structures. 


Consider a truss structure composed of a large number (N) of members of similar lengths, 
and where the shape distortion is due mostly to length errors (typically manufacturing 
errors). We assume that the length errors are uncorrelatcd and have the same standard 
deviation a c so that the covariance matrix E is 

S = a] I (28) 

To obtain the modes <j>j we expand the matrix B to the full set of finite-element degrees of 
freedom by adding rows and columns of zeroes denote it as B , treat it as a mass matrix and 
solve for the vibration modes from 

{K-u,]B)ij=0 (29) 

where 4>j denotes the j th B- vibration mode (excluding rigid body modes). The sensed 

components of <f>j are extracted as 

4>j - T<t> } (30) 

where T is a Boolean matrix. Note that 

B = T t BT (31) 

Because the rank of B is smaller than or equal to n 3 , only n s or less of the frequencies u>j 
have finite values. The infinite frequencies correspond to zero-rms modes satisfying 

4> T B<j>j = 4> T jB<t>j = 0 (32) 

To calculate the distortion associated with length error A/, in the ith member (corresponding 
to e, = 1) we apply a force vector Tq to the truss and solve 

/y'V\ = A/./q (33) 

where /q consists of a pair of forces colinear with the ith member of magnitude 

/q = {EA)i/li (34) 
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where ( EA)i and /, denote the axial rigidity and length, respectively, of the eth member. We 
can then extract the sensed components of 0,- 

0, = T 0, (35) 

The displacement vector can be written in terms of modal coordinates 

0, = (36) 

j=i 

where 

0 lj = 0 J T 2?0, = 0j50 t (37) 

Equation (37) implies that the modal coordinates of 0; and of 0,- are the same. This is a 
result of the fact that the modes beyond n, do not contribute to the rms 

n. 

(0, 2 )rm, = 0, T i?0, = 0, T W = 0* (38) 

>= 1 

The modal coordinates can be calculated without solving first Eq. (33) by premultiplying it 
by 0J and using Eq. (36) to obtain 

0jtf^0.,0, = A/.0JF, (39) 

<=i 

But 

%K4t = u o*6 jt (40) 

where is the Kronecker delta. Altogether we get 

0 U = A/,0j/'Vs 2 OH) 

We can simplify expressions by rewriting Eq. (40) as 

N 

<i>jK<i>( = Krft (42) 

1 = 1 

where /\, is the stifFness matrix of the ith element. Then 

N 

VjK~4>i = ^(0jF t )(0[F,)/,/(£;A) < = (43) 


8 



Finally using Eqs. (13), (28) and (41) we get 

N 2 N 

(C**),, = V = ,-fe E A, ?WJ S) (44) 

' UJrLUf . 

«=1 •> * 1=1 

Using Eqs. (14), (19), (25) and (44) we can calculate the optimum distortion ration ( g 2 ) 
opt. However, further simplification is possible since in comparing Eqs. (43) and (44) we 
see that if 

A if(EA)iHi = c (45) 


then 

(<?«),< = (46) 

U t 

This happens, in particular, when all the members of the truss have identical nominal prop- 
erties. In this case since C ^ is diagonal its eigenvalues are equal to the diagonal elements 
and Eq. (25) becomes 

= “4 £ "2 ( 47 ) 

jrrn+l i=n+ 1 i 

and the distortion ratio of Eq. (12) becomes 


(g 2 ) oP t = 


n> i 

j=n+ 1 J 
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E 4 ? 

j=i > 


(48) 


Note that in this case the eigenvectors of can be taken to be the unit coordinate vec- 
tors. That is, for the case of a truss satisfying Eq. (45) the optimum actuators produce a 
deformation state which is a linear combination of the first n B- vibration modes. 


Application to 55m Antenna Truss. 

The truss support structure for a 55-m radiometer antenna shown in Figure 1 was used 
as an example. The reflector is built up from tetrahedral truss modules, and consists of 420 
members connected at 109 joints. The upper surface and the r.m.s. of the vertical motion of 
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this surface was used as a measure of the error. This corresponds to the vector ip consisting 
of the vertical displacements of the 61 joints on the upper surface, and the matrix B equal 
to ( 1 /x/6l)I. 

The truss elements of the upper and lower surface have similar lengths of about 315 in. 
with nominal length variations of about one percent. The diagonal elements connecting the 
lower and upper surface are about 180 inches long with nominal length variations of about 
3-4 percent. This antenna structure was analyzed in Ref. 3 assuming that shape errors are 
introduced due to manufacturing length errors from the nominal values. The length errors 
were assumed to be uncorrelated and of equal standard deviation a c corresponding to all A/, 
being equal, and to E given by Eq. (28). The actuators were assumed to be mechanisms 
(such as screws) embedded in some elements which can be used to change the length of these 
elements. 

Two methods were used to obtain optimum actuator positions on the upper surface of the 
antenna in Ref. 3. First an approximate procedure that permitted actuators to be located 
anywhere in space(that is not restricted to truss member location) was used. This procedure 
employed the conjugate gradient method to minimize E(6f ms ). The second procedure was 
heuristic integer programming approach called ESPS (Exhaustive Single- Point-Substitution, 
Ref. 4). 

The optimum actuator locations for the cases of three and six actuators, are shown in 
Figures 2 and 3, respectively. The corresponding optimal values of the distortion ratio g, 
Eq. (12), are compared in Table 1 with the lower bound on g predicted by Eq. (48). It is 
seen that the approximate optimum (which permits actuators even where there are no truss 
elements) comes close to the lower bound, with the ESPS optimum not too far behind. This 
good agreement is achieved in spite of the fact that the assumption of Eq.(45) is not met 
exactly (with small variations for surface elements and a large discrepancy for the diagonal 
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elements). 


A clieck on the effectiveness index of the individual actuators used in the ESPS solution 
was performed next. For the three actuator case < 7 ,- ranged between 0.130 and 0.974 for the 
420 elements of the truss. The elements selected by the ESPS procedure as actuators, on 
the other hand, had effectiveness indices 9 , beween 0.961 and 0.974, indicating that they 
produce a displacement field overwhelmingly represented by the highest 3 modes of C It 
is interesting to note that only 29 of the 420 members had < 7 , larger than 0.95 so that the 
selection of actuator locations could have been greatly simplified. For the six actuator case 
(ji ranged between 0.164 and 0.980 with the ESPS procedure selecting actuators in the range 
0.966 to 0.980. Of the 420 members 36 had < 7 ,- larger than 0.960 and 45 had < 7 ,- larger than 
0.95. It should be noted that it is not reasonable to simply select the actuators with the 
highest effectiveness indices. It is important that the actuators produce displacement field 
that are not linearly dependent or close to being linearly dependent (see Ref. 8 ). 

These results have three implications. First they validate the procedures used in Ref. 
3 for obtaining good actuator positions. Second they indicate that Eq. (48) can be used 
to assess the degree of shape control possible with a given structure and a given number of 
actuators. Therefore, it can be employed in a structural redesign of the truss to improve 
shape controllability. Third, they show that the effectiveness index of Eq. (27) can be used 
to screen actuator locations and greatly reduce the number of actuator locations that need 
to be considered in the search for the optimum locations. For example, for the six actuator 
case, if we consider only the 45 members with < 7 , greater than 0.95, the number of possible 
configurations drops from 7 x 10 12 to 12 x 10 6 . 

Concluding Remarks. 

This work investigated the correction of shape distortion due to errors n structural sizes 
which are zero-mean normally distributed random variables. A bound on the maximum 
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improvement in the expected value of the root- mean-square shape error was obtained as well 
an effectiveness index of a given actuator. The results were specialized to a simple form 
for truss structures composed of nominally identical members. The bound was tested on 
a 55-m radiometer antenna truss structure. It was found that previously obtained results 
for optimum actuators had a performance close to the bound obtained here. These results 
raise the possibility that the bound obtained here will be useful in antenna truss design for 
improved shape controllability. Furthermore, it was verified that the optimum actuators 
found previously rank very high in terms of the effectiveness index developed here. The 
effectiveness index can therefore be used to screen potential actuator locations and reduce 
drastically the number of truss members that need to be considered in the selection for 
actuator placement. 
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Tabic 1 


Comparison of lower bound of distortion ratio with values achieved by optimization of 
actuator locations 


Number of 

lower bound 

ESPS 

approximate 

actuators 

Eq. (48) 

optimum 

optimum (Ref.3) 

3 

0.6384 

0.6919 

0.6595 

6 

0.5346 

0.5967 

0.5510 

9 

0.4617 

0.5271 


12 

0.4049 

0.4710 


15 

0.3531 

0,1061 


18 

0.3044 

0.3574 




Figure 1: Side and Top Views of Tetrahedral Truss 

Antenna Reflector 
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Figure 2: Optimum Upper Surface Actuator Configurations - 3 Actuators 



Continuous optimum, g= 0.6595 

hmi Best near feasible design, 
g = 0.7028 

ESPS Solution g = 0.6919 
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Figure 3: Optimum Upper Surface Actuator Configurations - 6 Actuators 



Continuous optimum 
g=0.55IO 

Best near by feasible design, 
g =0.6039 


ESPS Solution g = 0.5967 
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